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ABSTRACT 



We present a new global model of the solar corona, including the low corona, 
the transition region and the top of chromosphere. The realistic 3D magnetic 
field is simulated using the data from the photospheric magnetic field measure- 
ments. The distinctive feature of the new model is incorporating the MHD 
Alfven wave turbulence. We assume this turbulence and its non-linear dissipa- 
tion to be the only momentum and energy source for heating the coronal plasma 
and driving the solar wind. The difference between the turbulence dissipation 
efficiency in coronal holes and that in closed field regions is because the non- 
linear cascade rate degrades in strongly anisotropic turbulence in coronal holes 
(no inward propagating wave), thus resulting in colder coronal holes with the 
bi-modal solar wind originating from them. The detailed presentation of the 
theoretical model is illustrated with the synthetic images for multi-wavelength 
EUV emission compared with the observations from SDO AIA and Stereo EUVI 
instruments for the Carrington rotation 2107. 

Subject headings: MHD — Alfven waves — turbulence — coronal heating — heating 
function 



Introduction 



Results from Hinode observations have recently upped the estimates for the MHD wave 



energetics in the solar corona (De Pontieu et al. 2007). Observed in the chromosphere 



the magnetic field perturbations appear to be so powerful that even a 10 -i- 20% fraction 
of them, which propagate out from the Sun, carry a large enough energy to heat the solar 
corona and accelerate the solar wind. Even more promising, with the launch of the Solar 
Dynamics Observatory, we are now beginning to see observational hints of these ubiquitous 
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waves in the transition region and low corona {Mcintosh et al. 2011) 



Even before these encouraging observations had been obtained, models which 
incorporated, or even were entirely based upon Alfven wave turbulence as the momentum 
and energy source were developed to describe the solar wind and coronal heating. 
Nowadays, developing the turbulence- driven global space weather model becomes a problem 
tempting to be solved. 



1.1. Solar wind: can the turbulence-driven model compete with the 

semi-empirical one? 



In Usmanov et al. (2000), a three-dimensional (3D) model for the solar wind was 



suggested in which the Alfven wave turbulence pressure served to accelerate the solar wind. 
The solar wind bi-modal structure as observed by Ulysses was successfully reproduced in 
the numerical simulation. However, the quantitative agreement of this model with long 
history of the solar wind observations at 1 AU is not good enough. The same criticism 
seems to be applicable to the more refined and physics-based Alfven-wave-driven models 



of the solar wind Suzuki and Inutsuka (2005); Verdini et al. (2010); Osman et al. (2011) 



Therefore, the semi-empirical approach so far is better suited for global space weather 



simulations. The most popular parameterization was adopted by Arge and Pizzo (2000) in 



their Wang-Sheeley-Arge (WSA) model, which well describes the solar wind parameters 
at 1 AU. The semi-empirical "synoptic" formulae for the solar wind speed employ the 
solar magnetogram and the properties of the magnetic lines of the potential magnetic 
field as recovered from the synoptic magnetogram data. The empirical dependence of the 
solar wind properties from two input parameters: 6', which is the angular distance from 
the solar wind origin point till the nearest coronal hole boundary, and /cxp, which is the 
expansion factor, /exp = |B^=/jg |i?^/[|B/{=2.5_RQ 1(2.57?^)^], is derived from the following 
two assumptions and observations. First, the slow solar wind is assumed to originate 
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from the coronal hole boundary (small values of 6), whereas the fast solar wind originates 
from the central part of the coronal hole (large values of 6). Second, small coronal holes 
(having large values of /) usually produce slower solar wind compared to large coronal 
holes (having small values of /). 



In Cohen et al. (2007), we coupled the semi-empirical WSA formulae to the global 



three-dimensional model for the solar corona and inner heliosphere within the Space 



Weather Modeling Framework (SWMF) ( Toth et al. 2005). The WSA formulae were used 



as the boundary condition for the model which had been coupled to the MHD simulator 
via the varied polytropic index distribution (see Roussev et al. ( 2003b[ )). However, the 
physics of the Alfven wave turbulence has almost no intersection with the semi-empirical 
model of WSA. 

In order to he competitive with semi- empirical model of the solar wind, the turbulence- 
driven model should quantitatively reproduce the solar wind variation from the coronal hole 
boundary to the coronal hole central part, similar to that parameterized by the expansion 
factor. 



1.2. Coronal heating: can the turbulence-based model compete with ad hoc 

heating functions? 



Proceeding from the solar wind to the coronal heating, we can see again the 
disconnection between the physics-based models for the turbulent heating in the corona. 



one of the most advanced models of this kind being recently described in Cranmer (2010), 



on one hand, and well established models with semi-empirical heating function, such as 



that presented in Lionello et al. (2001); Riley et al. (2006); Titov et al. (2008); Lionello 



et al. (2009); Downs et al. (2010), on the other hand. The heating function is applied to 



power the plasma in the solar corona with some heating rate, the distribution of which is 
given by some functions chosen not from some deep physical considerations, but chosen in 
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an ad-hoc manner to provide a better agreement with observations. 

These heating functions should be apphed in conjunction with thermodynamic 
energy equation(s) in the Lower Corona (LC) model. A direct accounting is needed 
of the transition region between the chromosphere and the corona, where non-MHD 
thermodynamic terms of energy transport, such as electron heat conduction, radiative 



losses, and coronal heating all become important. In Downs et al. (2010) we used this 
thermodynamic MHD model to explore empirical parameterizations of coronal heating 
in the context of realistic 3D magnetic structures observed in the EUV and soft-X-Rays 
on Aug 27, 1996. Through direct comparison of synthetic observables to observations we 
demonstrated that this model can effectively capture the interplay between coronal heating 
and electron heat conduction. 

The plasma parameters distribution obtained in this way may be successfully applied 
to generate synthetic EUV and X-ray images, which appear to be in a good agreement 
with those obtained with the KIT telescope onboard SOHO and SXT onboard Yohkoh 
(up to 2001). The observation synthesis capability has been extended to the major low 
corona imaging instruments available in space, namely STEREO/EUVI, SDO/AIA, and 
Hinode/XRT. The best agreement with observations is achieved with ad hoc heating 
functions (such as the unsigned-flux-based heating model as discussed in Section 2.1). In 
order to be competitive with ad hoc coronal heating models, the turbulence- driven model 
should quantitatively reproduce some successful heating function, that is the realistic wave 
dissipation should provide the same heating rate, as, for example, the unsigned-fiux-based 



heating mode/ Abbett {2001). 
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1.3. Towards a Global Alfven- Wave-Turbulence-Driven MHD Model of the 

Solar Corona and Solar Wind 

In addition to tlie mentioned disconnection between pliysics-based and observations- 
driven models, there are two more contradictions to comment on. First, the quantitative 
Alfven-wave-driven models for the solar wind and for the coronal heating do not well 
conform with each other. For example, the Cranmer's model of the coronal heating 



{Cranmer 2010), once applied to realistic 3D global model does not give realistic plasma 
parameters in the solar wind region (in the coronal holes). Second, the increased estimates 
for the Alfven wave turbulence energetics in the Sun's proximity require us to revisit the 
models for the evolution of turbulence while the solar wind propagates towards 1 AU. 
The revisited model should account for both the observed level and frequency spectrum 
of turbulence at 1 AU and the solar wind ion temperature resulting from the turbulence 
dissipation. 

With the advent of modern computational tools it is now becoming the norm to 
employ detailed 3D computer models as empirical tools that directly account for the 
inhomogeneous nature of the Sun-Heliosphere environment. The key advantage of this 
approach is the ability to compare and validate model results through direct comparisons 
to all kinds of the observational data listed above: for the solar wind and turbulence 
characteristics at 1 AU and for EUV and X-ray images of the solar corona. 



1.4. Goal and Content of the Paper 

In the current paper we get rid of the ad hoc heating functions and parameterize the 
coronal heating in the LC in terms of the Alfven wave turbulence dissipation. 

The theoretical model for this approach is summarized in Section 2. In subsections 2.1, 
2.3 we demonstrate the possibility to parameterize a popular and successful semi-empirical 
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heating model based on unsigned magnetic flux in terms of the Alfven wave turbulence 
dissipation. While choosing the way to parameterize heating, one may try to vary both 
the boundary condition for the wave energy flux (the Poynting flux) inflowing to the solar 
corona from the solar surface and the dissipation length for turbulence. In subsection 
2.2 we show that the boundary condition is strongly restricted which reduces the model 
uncertainty. With the pre-specified boundary condition, the drastic difference in the plasma 
heating mechanism between closed field region and coronal holes should be caused by a 
difference in the wave dissipation efficiency. In subsection 2.4 we suggest and describe the 
physical mechanism of a nonlinear interaction between oppositely propagating waves. The 
degraded intensity of the inward propagating waves may be responsible for the reduction 
in the turbulence dissipation rate in the coronal holes thus resulting in the bimodal solar 
wind structure. 

In Section 3 we summarize a computational model and the code we use to simulate 
the state of the solar corona. Section 4 presents the simulation results for CR2107 and 
their comparison with EUV images. In the Conclusion we discuss the plans for a future. 



2. Coronal Heating and Its Parameterization via the Alfven Wave Turbulence 



2.1. Coronal Heating Model Based on the Unsigned Magnetic Flux 



Among the possible choices for the volumetric heating function examined in Downs 



et al. (2010), the most advantageous one reproduces both EUV and soft X-ray observations 



was that adopted earlier on by Abbett (2007). This heating function is based on the scaling 



law obtained by Pevtsov et al. (2003). This law establishes the power-law relationship 



between the total heating power, E = J edV, integrated over the plasma volume above 
IRq, and the unsigned open magnetic flux, $ = J \^Rq ■ dS|, integrated over the 
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photospheric surface, i.e., 

erf\/ = const $, $ = f |B-dS|. (1) 

Here R is the heliocentric distance and Rq is the solar radius. In the original work by 
Pevtsov et al. (2003), the relationship was found not to be linear: Ex ~ 0.894$^'^^^^ 



[CGSE]. Note, however, that the observational data used to derive the scaling law were the 
X-Ray total luminosity. Ex-, which is only a small fraction of the total heating power, E. 
There are also losses due to the electron heat conduction, radiation in spectral ranges other 
than X-rays, and also solar wind expansion. As a result, a scaling factor that relates Ex 
to E is required, C = Ex/ E ^(1^2)- 10^^. If one assumes that Ex/ E ratio is somewhat 
elevated during solar maximum, or near active regions where the average magnetic field 
strength, |_B|, is higher, so that C oc ^o-^^ss^ then we can state that Eq.([T]) is equivalent 



to that presented in Abbett (2007); Pevtsov et al. (2003). On one hand, this relationship 



quantifies the well-known fact that during solar maximum, or near active regions, the 
average B is higher. On the other hand, the coronal heating is more intense in this case. 

The constraint given by Eq.([T]) is not sufficient to establish the 3D distribution of 
the heating function. In Abbeti\ ( |2007[ ), the heating function was scaled with the magnetic 



field: e oc |B|, with the constant factor being chosen to satisfy a power-law scaling 
similar to Eq.([T]). In \Downs et al. (2010) an exponential envelope has been adopted: 



e oc |B|exp[— (i? — Rq)/L], where the dissipation length was chosen L 40 Mm. The 
latter formula, however, appears to be applicable only for closed field regions away from 
active regions. In coronal holes, in turn, the dissipation length was chosen to be by a factor 
of ten greater. In all the cases, the sophisticated envelope function had to be integrated 
over the volume in order to obtain the common constant factor satisfying Eq.([T|. 

Below we demonstrate that a model of coronal heating that both satisfies Eq.([T]) and 
provides the spatial distribution of the heating envelope similar to that discussed above 
can be realized assuming Alfven waves being the main heating source the solar corona. 
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Without going into details about the wave absorption and non-linear conversion and 
transport, we can construct the heating function by: (i) imposing a boundary condition 
for the Poynting flux of Alfven wave energy; and, (ii) adopting an absorption mechanisms 
for the waves that would result in the desired distribution of energy deposition from waves 
into the coronal plasma. 



2.2. "Percolation" Boundary Condition for the Poynting Flux at the 

Photosphere 

Assuming that the Alfven wave turbulence is the main source for heating the coronal 
plasma, here we discuss a choice of boundary condition for their Poynting vector, P at the 
solar surface. We consider that the local value of the Pointing flux at the photosphere, 
Prq, scales with the magnetic field magnitude, iBl^g as: 

= const |B|rq. (2) 

With this choice for P, Eq.([T]) is automatically satisfied if one assumes that the entire energy 
of the Alfven waves traveling from the photosphere to the solar corona: Jj^^j^^ l^j^l-^ ' ^'^1' 
is deposited into the coronal plasma, /^>^g e dV. Then: 

E = I edV= I ^^|B-dS| = 



R>Rq JR=Rq \^\rq 



|B|rq Jr=Rq |B|ft 



meaning that Eq.Q is fulfilled, and that the exact same constant is present in Eqs.( Ip ): 

E Prp, 



$ |B| 



Rp 



(3) 



However, the relationship given by Eq.Q is not the only way to arrive at Eq.g. It 
can be also thought of as the percolation boundary condition for the Poynting flux of the 
waves. Assume that all the magnetic field in the solar corona and all the Alfven wave 
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turbulence percolate to the solar corona from a depth at the solar interior where both the 
magnetic field intensity (presumably 50 — 70 kG) and the Alfven wave Poynting fiux are 
very high compared to their photospheric values (except in sunspots). In this case, both 
the magnetic field fiux and the wave energy fiux "propagate" along the fiux tube, and they 
both vary inversely proportional to the fiux tube cross-section. (Note that the growth of 
the fiux tube cross-section from the solar interior towards the solar surface may be quite 
large and is hardly predictable.) As a consequence, the local value of the magnetic field 
at the photosphere is orders of magnitude smaller than that in the solar interior, and its 
variation along the solar surface may be quite large. The same behavior should be true for 
the Poynting fiux. However, the ratio of the Poynting fiux to the magnetic field intensity 
remains constant along the magnetic tube, and that is exactly what is assumed in Eq.(|2]). 

Surprisingly, in addition to these two considerations, we can cite four different groups 
of authors, which arrive at the same condition as in Eq.([2]) by quite different reasons (and 
stemming from quite different observations). First, the assumption as in Eq.([2]) applied to 
the radial components of the Poynting flux, Pr, and the magnetic fleld, Br, is a keystone 
of the Fisk theory for the solar wind (see \Fis^ ( |1996D ; \Fisk et aL\ ( |1999b|[a| ; \Fis^ ( |2001[ ); 
Fisk and Schwadron (2001). Then, \Farrugia et al. (1997) (and the works cited therein) 
formulated an "abnormal" adiabatic expansion law relating the (turbulent) energy density, 
w to the mass density, p: w oc p^^^ for the coronal matter, based on numerous observations 
relating to the coronal plasma and the CME ejecta. The suggested adiabatic law has an 
abnormal value of the adiabatic index: w oc p'^ , '-f = ^ < 1. Here, assuming a plasma 
dominated by turbulence, p is the mass density and w could be the energy density of the 
Alfven waves {w ~ {5BY, with 5B being the irregular magnetic fleld). As long as for the 
Alfven waves P = Vaw, and the Alfven wave speed scales with the density and magnetic 



fleld as Va = 
adiabatic law, 



|B| 

/mop 



oc Bp ' , Eq.(2) may be also rewritten in the form of the "abnormal" 



Wr^ = const 



P 
IBI 



(4) 



As long as media with adiabatic index less than one are not thermodynamically stable, it 
is hardly instructive to interpret Eq.([2]) in this manner. However, this unexpected support 
for Eq.(|2| is worth mentioning here. 

Another example supporting Eq.(|2| is the paper by Suzuki (2006), which stems from 
the Wang-Sheeley-Arge model, and proposes a semi-empirical quantitative model for the 
solar wind. In their work, it is assumed that the solar wind is mostly powered by the 
Alfven wave turbulence, and the input parameter in the model is the constant value of the 
following average: < 5B_|_ ■ 5v_|_ >^ 0.83 Tm/s at the solar surface. Here, 5B_|_ and 5v_|_ are 
the turbulent magnetic field and the turbulent velocity pulsation, respectively. These are 
both being orthogonal to the regular magnetic field. It is straightforward to demonstrate 
that the above assumption is equivalent to Eq.(|2]): 



IBI w 



|B 






B| 



Bl 



From here, one can see that the Alfven wave turbulence may be employed to reproduce 
the observed solar wind parameters with the use of a boundary condition similar to 
that given by Eq.([2]). In this case, the constraint for the constant factor in Eqs.(T]|2) is 
(P/|B|)i?Q > (0.83 Tm/s) //^o ^ 7 ■ 10^ W/(m2T). This is the lower bound as long as in the 



model of Suzuki (2006) there is no energy loss mechanism from the solar wind plasma (e.g. 
radiation, heat conduction, etc.). Therefore, the realistic heating rate should be higher in 
order to balance these losses. Nevertheless, the estimate is very close to that which follows 



from Pevtsov et al. (2003); Abbett (2007); Downs et al. (2010). Specifically, for the choice 



C = 2 ■ 10-2 (|B|[G])"-'^'', where the average magnetic field intensity over the solar surface 



is introduced, the estimate of the Pointing-flux-to-field ratio is as follows: 



E Ex 



0.894$^-i488 



2 ■ 10-2|B|0-i488<|, 
or, in all equivalent forms: 



.2 -jO.1488 



[CGSE] ) 



1.1 • 10 



W 



m^T ' 
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^ 6.6 ■ 10~^W/m^ ^ 
Va V 3 ■ 10-iikg/m' 



w = — ^ 6.6- 10~^W/m%/ ^ _ ^ 3 , < (5Bx ■ (5vx >^ 1.4Tm/s (5) 



< (5v_L ■ S\± >^ (15km/s 



,2. 



'3 ■ 10-iikg/m 



3 



where we chose the value of mass density 3 ■ lO^^^kg/m^ (A^^ ^ 2 • lO^^m^^) corresponding 
to the top of the chromosphere (near the transition region), to compare with the Hinode 
observations of the turbulent velocities ~ 15 km/s. Thus, a variety of models, which are 
validated against EUV observations. X-ray radiance measurements, as well as in-situ solar 
wind observations are all in agreement with Eq.(|2]), with a rather narrow range for the 
constant factor in the expression. 

Next, let us compare this boundary condition with a coronal heating model based 



on Alfven wave turbulence, such as the advanced model of Cranmer (2010) (see also 
the works cited therein). By comparing Eq.([2| with Eq.(29) from Cranmer (2010), one 
can see that Eq.(|2]) is adopted in their model not just as a boundary condition at the 
solar surface, but as a fundamental relationship valid throughout the solar corona (or, 
at least, for heliocentric distances I.IRq < R < 1.5Rq). The observational constraints 
for the constant factor in Eq.(|2| are somewhat weaker than the values discussed here 
(P/|B| ~ (0.5^0.8^ 1.0) ■ 10^W/(m^T)). However, the deviation is very small. In addition, 
there is some uncertainty in these observations (e.g., the Pointing flux is not measured, 
but rather the oscillating velocity). Also, the solar observations at I.I-R0 cannot be easily 
converted into a boundary condition at the photosphere. Furthermore, the dissipation 



length of 40 Mm assumed by Downs et al.\ (2010) is shorter than the difference between 



1.1 Rq (where the oscillations are observed) and 1 -R0 (where the boundary condition is 
imposed). As a consequence, it is likely that the turbulence observed at 1.1 Rq is already 
weakened by the wave absorption. 
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2.3. Heating Function Parameterized in Terms of Alfven Wave Dissipation 

In order to proceed to the 3D distribution of the heating function, we represent the 
volume integral of the heating function as the total of integrals over the flux tubes: 



E = j (^j edSii)^ di, 



with di being the length differential along the flux tube, and dS{i) being the expanding 
cross-section of the flux tube. Along the flux tube, the magnetic flux is conserved, and 
hence: dS{i)\B\{£) = const = |(B ■ dS)|. Again, assuming that the Alfven waves propagate 
only along the flux tube, their energy conservation is given by the following relationship: 

P{e + di)dS{e + di) - P{i)dS{i) = -edidS{i), hence: 

|B|(C) dt \\Bm) > 

By introducing the absorption length, L, such that e = Vaw/L = P/L, one can integrate 
Eq.([6]) along the magnetic fleld line: 



lBl^i?=RQ ^ \ Jo L 



Depending on the choice of the absorption length, one can obtain the following spatial 
distribution of the heating function: (i) for large L, the heating function is the same as 
that as used in Abbett ( 2007[ ), i.e., e oc |B|; (ii) for small L, the heating function decays 



exponentially, i.e., e oc exp(— J di/L). We now arrive at an important conclusion, that 
is, the most common heating functions may be parameterized in terms of an effective 
absorption coefficients for the Alfven wave turbulence! The desired absorption coeflicient, 
Va/L, may be directly included to the WKB equation for total wave intensities, w±, 
propagating parallel and antiparallel to the mean magnetic fleld vector, b = B/|B|, 



correspondingly ( Usmanov et al. 2000): 



^ + V ■ (u^± ± hVAW^) + \w^{W ■ u) = (8) 
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However, Eq.(|8|) can only be applied as a rather crude approximation. In the more 



refined model as derived in Sokolov et al. (2009)) on the left hand side of Eq.rtSl) the term 



proportional to the velocity divergence and spectral evolution in the convergent /divergent 
fiows should be properly accounted for. Even more important is the role played by the 
frequency dependence on the right hand side of Eq. (Isj) : 



^ + V ■ (uJ± ± hVAl±) - ||^(V ■ u) = -r(a;)J±, (9) 

where the equations are formulated for the wave energy density, I±, related to the unit 
volume and to the interval of the wave circular frequency, du, in the co- moving frame of 
reference, so that w± = I^du. Not only is the dependence of the dissipation coefficient, 
r(ci;), on the wave frequency significant (practically so that only for the highest wave 
frequencies is the absorption non-negligible), but also the nonlinear wave- wave interaction 
may occur in the form of the wave package upshift in frequency. As a result, the wave 
energy from low-scale perturbations fiows through the spectrum towards the short-scale 
end and once transferred to the shortest possible scales, the energy dissipates. 



2.3.1. Isotropic Kolmogorov Turbulence and Dissipation 

We found above that the unsigned-fiux-based heating function with the exponential 
envelope function can be realized as the dissipation in the Alfven wave turbulence. In 



{Downs et al. 2010) we saw that the unsigned-fiux-based heating function well reproduces 
EUV and X-ray images of the closed field regions (about the coronal holes, see below) 
with the choice for the dissipation length to be ~ 40 Mm. Here we discuss a choice for the 
dissipation length, L, based on the phenomenological turbulence theory, desiring that in 
the closed field region L would be about 40 Mm and do not vary too strongly, at least at 
low altitudes < 0.1i?Q ^ 70 Mm. 

Consider extra terms in the governing equations for the Alfven, turbulent energy 
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density, accounting for the wave-wave 'cascade' rate. Only the cascade-describing terms 
are kept. For the closed field region we employ the key assumption about the isotropic 
turbulence: 

/+ = /_ = /. (10) 

For simplicity, assume that T^wlu) characterizes the rate of conversion of two wave 
"photons" of the frequency of ~ w to a single photon with the frequency, ~ 2u. As 
long as the wave energy conserves in this processes, this may be thought of as the 
conservative advection in the frequency space. The positive speed of this advection being 
approximately ~ uT^wi^), as long as with the rate, F^^^, the wave frequency is upshifted 
by [(~ 2u) — (~ u)]. So, we have a phenomenological equation as follows: 

^ + -^(uTUu)I) = -T{co)I (11) 

The dimensionless ratio, T^^/u, may be parameterized in terms of another dimensionless 
ratio, uI/B^, B being the magnetic field magnitude. For two most natural estimates. 



F^^ oc uj ■ {ujI /B"^) oc ijj{6B/By or F^^,^ oc uj^/ujljB'^ oc uj{5B/B), we arrive at the 



Kraichnan spectrum or the Kolmogorov spectrum as the steady-state solution of Eq.(ll) 



throughout the inertial range of frequencies (we neglect wave absorption in this range): 





d_ 
did 



for / oc cu 3 ' 2 . We chose the assumption of the Kolmogorov spectrum, so that 

r ^ M_ 

J- WW T r \ 

VaV p 



On integrating Eq.(ll ) over frequency from some frequency value within the inertial range, 
cOin, to infinity with neglecting wave dissipation within the inertial range, we find: 

n.)M.J^M) . (12) 

The equation is valid for any choice of Uin within the inertial range, as long as the right 
hand side is constant for / oc u}~^^^ and F(co'j„) = 0. The energy of the Kolmogorov 
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spectrum, (x. J co ^^^dco, is concentrated at low frequency, therefore, at the least possible 
~ Umin oue cau estimate {u}I)uj=urm„ ~ and 



/ L ioojlauj (X — — A —w, 

J^ir. Va y p 



Following Hollweg (1986), for an isotropic Alfven wave turbulence we accept 



1 [ojjldoj ~ —w = w, [Id) 

so that the total volumetric heating rate due to the turbulent energy dissipation equals: 



1 Iw = = w+ . , ,^ ^, 

— {W-+W+). (14) 



The introduction of Lj_ oc l/kj_, which is the transverse correlation length of turbulence 
instead of u/Va = k\\, which is the wave vector for the Alfven wave propagating along 
the magnetic field line, makes the cascade theory better suited to the nature of purely 
transverse Alfven waves. For such waves the ratio of the transverse turbulent velocity to the 
longitudinal wavelength could hardly characterize the rate of nonlinear wave dissipation. 
In more refined turbulent theory (which we do not apply and even do not review here) the 
spectral energy should be introduced separately for parallel and transverse wave vectors. 
However, within our phenomenological approach we assume that kj_ oc k\\ and combine, 
once needed, the spectral energy distribution over u = VAk\\ with the transverse corelation 
length, to quantify the non-linear dissipation rate. 

To close the model and to compare its predictions with the unsigned-flux-based heating 



function, we chose, again following Hollweg (1986), the scaling law for the transverse 
correlation length: 

,^,_i/2 , 75km ■Ti/2 
L_L oc |B| L_L = ^= — . (15) 

In the present work we use the following range of values for the empirical constant: 
L±^/\B\- 

50 km ■ T^/2 < L , ^AB\ < 100 km • T^/^ 
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Hollweg (1986). The 



in the numerical estimates below we use L± ^/\B\ ~ 75 km ■ T^/^ ^^^^ 
scaling law |B|L5^ = const is well compatible with the "percolation hypothesis" noticed 
above: while waves propagating from the solar interior along the flux tube towards the solar 
surface it is natural to assume that varies proportionally to the flux tube cross-section 

To discuss the value and the variability of the spatial scale in the coronal heating, if 
the latter is related to the Alfven wave turbulence dissipation, compare Eqs.( l4p5 ) with 



the representation of the wave dissipation as in the right hand side of Eq.(|8j). From this 
comparison we find: Va/L = y^w/ p/L±. At the corona base we can employ Eq.(j4]) to 
evaluate the boundary values of w and the estimate for L in this region reads: 



L = , ^ ^ J ^^40 Mm. 



We see that with the choice of the above estimates for the transverse correlation length 
and for the boundary condition for the Poynting flux the estimate of the dissipation length 
at the coronal base has only weak scaling (as a square root) with the Alfven speed. The 
latter does not vary too much at the coronal base, as long as Va oc |B|/yp and the density 
correlates with the magnetic field intensity, being larger in the active regions with strong 
magnetic fields. We also see that for a reasonable estimate for the Alfven speed at the 
coronal base, Va ~ (200 — 500) km/s, the dissipation length is close to the ad hoc value 



as used in Downs et al. (2010). The dissipation length is not used below directly, as 
long as Eqs.( l4p5 ) do not include it. However, the capability to reproduce the popular 
unsigned-fiux-based heating model with the popular envelop function is important. We can 
now strengthen the conclusion above in the following manner: the choice of a boundary 
condition for the Poynting flux together with the estimates of the non-linear dissipation in 
an isotropic Kolmogorov turbulence allows us the reproduce in detail the unsigned- flux-based 
ad hoc heating model for the closed field region. 
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2.4. Coronal Holes and the Solar Wind Model 

An important ad hoc approach while applying the heating function is the capability 
to use strongly different functions in the coronal holes and in the closed field regions. 



In Downs et al. (2010); Lionello et al. (2009) in the coronal holes the following heating 



function was applied: 

e = 5 ■ 10"^ exp[-(i? - Rq) / {<d.7RQ)\ erg/cm^s. 

Comparing with the unsigned-fiux-based model for the closed field region, the spatial scale 
of the heating function in the coronal holes is significantly longer than ~ 40 Mm which we 
used above. 

Discuss the possibility to implement such heating model via the Alfven wave 
absorption. First, we do not assume that in the coronal holes the Poynting flux at the 
corona base is different from that in the closed fleld region and attribute the drastic 
difference in the heating functions to the difference in the dissipation rate only. Second, 
comparing the integral of the heating function over the coronal hole volume with the 
accepted value of the Poynting flux for the turbulence, one can flnd that only a few percents 
of the Poynting flux is absorbed in the coronal holes within the heliocentric distance range. 
Rq<R< 2.5Rq. 

Comparing with the closed field region, where almost all wave energy is absorbed, one 
can conclude that the wave energy dissipation rate in coronal holes should include a small 
factor of the order of a few percents compared with Eqs. [l^fT^. 



2.4-1 ■ Anisotropic Kolmogorov Turbulence and Dissipation 

The drastic decrease in the wave dissipation rate in the coronal holes comes naturally 
if we take into account a strong turbulence anisotropy in the coronal holes: the outward 
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propagating waves should be much more intense than the waves propagating toward the 
Sun. 

Now we consider an anisotropic turbulence, such that /+ 7^ As long as the real 
cascade physics requires the presence of oppositely propagating waves, we require that the 
wave-wave interaction rate should be the function of minimum of w_ , in a strongly 
anisotropic turbulence. To achieve this, in the expression for the dissipation rate we employ 
a function, 2w^w+/ {w^ + w+) which comes to w = w_ = w^, for w_ = w^, but which 
tends to 2 min(iy_, at m.m{w^,w+) ^ max(w_, so that 



|B| / 2w_w+ 

[W- + w+j 



For the closed field region this expression is more realistic than that we used before, as 
long it correctly captures more intense heating in shorter loops, in which the intensities of 
the oppositely propagating waves are more uniform, hence, more isotropic, as compared 
to longer loops, in which the wave energy for counter-propagating waves may be strongly 
non-uniform and non-isotropic. However, for coronal holes this approximation is not 



sufficient as long as the WKB approximation, under which Eq.(ll) was derived is not 
accurate enough. Within the model we use here, in the coronal holes only outward 
propagating waves can exist, while in reality the wave outward propagation is accompanied 
with some reflection, resulting in appearance of waves propagating towards the Sun. Small 
ratio of the reflected wave amplitude to the bulk wave amplitude is described by the 
reflection coefficient, Cred- Therefore, the ratio min((w_, w+)/max(w_,w+) should be 
floored with C^^g^, giving an ultimate expression for the turbulent dissipation rate we use in 
the present work: 



IBI 



,max(w_,w+) / 2w-W^ 

L'refl \ I 1- 



p y {w_ + w+) p 



{w^+w+) (16) 



- 20 - 



2.4-2. Estimates for the Reflection Coefficient 

In the present research we focus on the study of the Lower Corona and do not compare 
the solar wind generation and its propagation towards 1 AU. Therefore, for our present 
purpose the estimate 

Crefl = const ~ 0.01 ^0.1 

is sufficient. 

Having in mind to develop this approach in the forthcoming publications, discuss 
briefly, how the idea of the WSA semi-empirical model can be implemented within the 
Alfven-wave-turbulence-driven model. As discussed above, for the WSA model the speed 
of the solar wind originating from the given magnetic field line footpoint may be found 
from two characteristics of the line - the expansion factor, /exp, and the distance to the 
coronal hole boundary, 9: 

In the original version of the WSA model, for example, the solar wind speed, was a function 
of the expansion factor only: — A/ f^^^, where the unknown constants, A and 5 were 
chosen for better fitting the observed solar wind speed at 1 AU. 

Once formulated in terms of the Alfven wave absorption, the model can no longer 
employ the existing WSA formulae for the solar wind, as long as there is no simple 
relationship between the wave absorption and the solar wind speed. However, we can keep 
using the idea of the WSA model and fit the observed solar wind parameters at 1 AU by 
properly choosing a formula for Crefl that depends on 9 and /exp- The first choice to be 
tested is: 

aefl = A . (17) 

where the unknown constants will be chosen to better fit the observations. 

Note that there is some reasoning in favor of this approach within the wave-based 
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model for the solar wind, as long as the expansion factor for a given magnetic field line 
is a good characteristic of the reflection coefficient for the waves, propagating along this 
line. Indeed, the expansion factor is larger for strongly bent magnetic field lines, but the 
reflection coefficients at such lines is also larger than on straight lines. At small distances 
from the coronal boundary the abrupt gradients of the plasma density also can be a reason 
for the increase in the reflection coefficient, resulting in the generally recognized opinion 
that the slow slow wind (higher wave absorption) has its origin near the coronal hole 
boundary (small 9). 

3. Computational Model 

In this section, we describe the governing equations to be solved numerically as well 
as the numerical tools we use in the numerical simulations. 

3.1. Governing Equations 

The model includes the standard MHD equations (non-specified denotations are as 
usually) : 



dp 
dt 



+ V-(pu)=0, 



(18) 



dB 

~dt 



+ V • (u (g) B - B (g) u) = 0, 



(19) 



^ + V.(pu^u-^)+v(p. + n + 



2/io 



+ 



2 




(20) 



with the full energy equations applied separately to ions 




(U ■ V)Pe + 



Nike 



(Te -Ti) + T {w. + w+) , 



(21) 
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and to electrons: 



+ V 



P. 



-u 



PeV ■ U = V (/tVTe) + 



f22) 



(9t ■ ■ V7-I 

In addition to the standard effects, the above equations account for a possible difference in 
the electron and ion temperatures, the electron heat conduction parallel to the magnetic 
field lines: 



K 



~B2 



-K.il, Kll OC 



rp5/2 



(23) 



the radiation energy loss from an optically thin plasma: 



Qrad = N,NiA{T,), 



(24) 



as well as the energy exchange between electron and ions parameterized via the relaxation 
time, Tei, as this is usually done. The Alfven wave turbulence pressure and dissipation 
rate is applied in the above equations. In particular simulations the heating due to the 



turbulence dissipation (see Eq 21) may be split between electron and ions and we can also 



use the total energy equation for electron and ions, to improve the computational efficiency. 
The turbulence propagation and dissipation are described within the WKB approximation: 



duj 1 

— h V ■ (uw± ± hVAW±) + -w±V ■ u = — rw±. 

ot 2 



(25) 



Summarize the above consideration and derivations regarding the dissipation rate, which 
we apply in the following form: 



|B| 



|B|) 



<-^rcfil / H 



2w_if_|_ 
(w- + w+) p 



(26) 



and with the following ranges for the parameters involved: 



[L^ ■ v1b|) = (50 ^ 100) km ■ T^/^ C.^a = 0.01 -^ 0.1, 



(27) 



We use the simplest equation of state for the coronal plasma with the polytropic index. 



7 



3' 
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The system of governing equations is solved numerically using BATS-R-US/SWMF 
code (see 3.3 below). The boundary condition for the Poynting flux (or for the intensity of 
the outgoing waves) is given by Eq.([5]). The boundary condition for the coronal magnetic 
field is taken from the full disk magnetogram. The boundary condition for the density and 
temperature requires more attention and is discussed in the following subsection. 



3.2. Chromosphere and transition region 

3.2.1. Chromosphere Boundary Condition 

Discuss the analytical solution of the the hydrostatic equilibrium and heat transfer 
equations with the exponential heating function, = Aexp{—x/L), The solution is as 
follows: 

ruigx 



Te = T, = To, Ne = Ni = Noexp 



kB{Te + Ti] 



Qh = Qrad = NlK{T,) = NlK{To) exp (^-^ j • (28) 

Here g = 274 m/s^ is the gravity acceleration near the solar surface, the direction of this 
acceleration being antiparallel to the x-axis and rrii is the proton mass. The two constants 
in the solution, Nq and Tq, which are the boundary values for the density and temperature 
correspondingly, are unambiguously related to the amplitude and decay length of the 
heating function: 

L='^^To-{30m/K), A = N^A{To). (29) 

Notice a very simple relationship for the exponential decay length for the heating function, 
which is the half of the barometric scale length for density: 2L = Lg = ksiTe + Ti)/{mig). 

The described solution satisfies the equation for the heat conduction as long as the 
heat transfer in the isothermic solution is absent and heating at each point exactly balances 
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the radiation cooling. The hydrostatic equihbrium is also maintained, as long as 




The suggested solution well describes the chromosphere. The short decay length of 
the heating function, which is equal to ~ 0.6 Mm for Tq = 2 ■ 10^ K may be, presumably, 
related to absorption of (magneto) acoustic turbulent waves, rapidly damping due to the 
wave- breaking effects. Physically, including this chromosphere heating function would 
imply that the temperature in chromosphere is elevated comparing to the photospheric 
temperatures due to some mechanism acting in the chromosphere itself. By no means this 
energy can be transported from the solar corona as long as the electron heat conduction 
rate at chromospheric temperatures is very low. 

As long as we do not apply such ad hoc short-scale heating function in the 
chromosphere, nor we include short-scale turbulent dissipation length, we set the boundary 

condition for the density and the temperature at the top of chromosphere: 

T,h = To = 2- 10" K, Nch = No^2- 10^^ m'^^ (30) 



3.2.2. Transition region 

The analytical solution for the transition region had been published many times. Here 
we focus on the following issues: (1) merging this solution to that for chromosphere; (2) the 
validity of zero-gravity approximation; and (3) the validity of the modified heat conduction 
model. 

The heat transfer equation for a steady state hydrogen plasma in a uniform magnetic 
field reads: 

£ ('^oT,^^'^) +Q^- ^"m) = 0- (31) 

Here Qh — r(w_ -|- w+) is the coronal heating function assumed to be constant at the 
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spatial scales typical for the transition region. Note that the coordinate is taken along the 
magnetic field line, not along the radial direction. 



On multiplying Eq.(31) by K,oTe^^{dTe/ds) and by integrating from the interface to 
chromosphere till a given point at a temperature, Tg, one can obtain: 

Here the product, NeTe, is assumed to be constant, therefore, it is separated from the 
integrand. For a given T^h the only parameter in the solution is {NfTf.)- In can be expressed 
at any point in terms of the local value of the heating fiux and the radiation loss integral: 



{NeTe 



7/2 
ch 



(33) 



The assumption of constant (NeTe) is fulfilled only if the effect of gravity is negligible. 
Quantitatively this condition is not trivial, as long as both the barometric scale and 
especially the heat conduction scale are the functions of temperature. The barometric 
scale may be approximated as Lg{Te) ~ Te ■ (60 m/K). The heat conduction scale, L^, can 
be estimated by noticing that within the large part of the transition region the radiation 
losses dominate over the heating function, therefore, they are balanced by heat conduction: 
KoTe^"^ ■ (Te/Ll) ~ Qr- Thus, the condition for neglecting gravity is: 



9/2 



MT.)«r..(60m/K)»L,.«,/^^^2|_. (34) 
In Fig. 1 we plot temperature dependencies, Lh{Te), Lg{Te), for (NeTe) = 10^° K/m^. We 



see that near the chromosphere boundary the approximation (34) works very well as long 
as the temperature changes with height are very abrupt. The increase in temperature till 
lO^K occurs at the length shorter than 0.1 Mm. This estimate agrees with the temperature 
profile recovered from observations of chromospheric lines (see, e.g., Fig.2 and Fig.8 in 



Avrett and Loeser (2008)). However, as the temperature further grows with the height. 
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the gravity effect on the tmeperature and density profiles becomes more significant. It 
becomes comparable with the heat conduction effect at ~ 4.5 ■ 10^ K, which can 
accepted as the coronal base temperature, so that the transition region corresponds to 
the temperature range from Tch ~ 2 • 10^ K till 4.5 ■ 10^ K, with a typical width being 
~ 10 Mm ^ Rs /70. The transition region solution merges to the chromoshere solution with 
no jump in pressure. The merging point in the chromosphere, therefore, is at the density of 
[NeTe) /Tch ~ 10^^ cm~^. The short heat conduction scale at the chromosphere temperature 
(see Figl) ensures that the heat flux from the solar corona across the transition region does 
not penetrate towards higher densities. 



It is known that Eq.(33) may be used to establish the boundary condition for density 



at the 'coronal base', however, we use another way to model the transition region by 



increasing artificially the heat conduction in the lower temperature range (see Abbett 



(2007)). Consider a following transformation of the temperature functions in Eqs.(31 32) 



1^0 fl^O, ds^ f ds, , r ^ r//, Qrad ^ Qrad/f, (35) 

with a common factor, / > 1. The equations does not change in this transformation and 
the only effect on the solution is that the temperature profile in the transition region 
becomes a factor of / wider. By applying the factor, / = (Tm/Tg)^/^ at Tch <Te< Tm, one 
can achieve that the heat conduction scale in this range is almost constant and is close to 
^ 1 Mm for a choice of T^ ^ 2.5 • 10^ K (see Fig. 1). 

It should be emphasized, however, that the choice of the temperature range to apply 



this transformation is highly restricted by the condition as in (34). In choosing a higher 
value of Tm the heat conduction scale at the chromospheric temperature exceeds the 
barometric scale in the chromosphere resulting in physically meaningless penetration of 
the coronal heat to the deeper chromosphere. The global model of the solar corona with 
this unphysical energy sink suffer from the reduced values of the coronal temperature and 
produces a visible distortion in the EUV and X-ray synthetic images. Thus, in formulating 



the transition region model we modify the heat conduction, the radiation loss rate and the 
wave dissipation rate and the maximal temperature for this modification does not exceed 
Tm ~ 2.5 ■ 10^ K. 



3.3. BATS-R-US and SWMF Codes 

The BATS-R-US (Block Adaptive Tree Solar Wind Roe-type Upwind Scheme) 
code has been developed at the UM. It solves the equations of ideal MHD — a system 
of eight equations describing the transport of mass, momentum, energy, and magnetic 



flux [Groth et al. 2000 Powell et al. 1999). This massively parallel code enables 



Sun-to-Earth simulations to be performed in near real-time when run on hundreds 



of processors on a supercomputer {Manchester et al. 2004b). The implementation of 



adaptive-mesh-refinement (AMR) in BATS-R-US allows orders of magnitude variation in 
numerical resolution within the computational domain. This is important for a global 
model of the solar magnetic fields in which one strives to resolve such structures as 
shocks, volumetric currents of flux ropes, electric current sheets in a 3D domain, which 
may extend out to hundreds of Rq. The use of AMR also enables us to resolve the fine 
structure of active regions on the Sun, which spawn CMEs. This is vitally necessary for 
incorporating high-resolution magnetic data into a numerical MHD model. In the context 
of solar-heliospheric physics, BATS-R-US has been utilized to model the global structure 



of the solar corona and solar wind [Roussev et al. 2003b; Cohen et al. 2007, 2008), as well 



as the initiation [Roussev et al. 2003a Jacobs et al. 2009) and propagation of idealized 



[Manchester et al. 2004a||b ) and not-so-idealized [Roussev et al. 2004, 2007) solar eruptions 



and associated SEP events [Sokolov et al. 2004) 



The SWMF (Space Weather Modeling Framework) is a high-performance computa- 
tional tool that has been developed at the University of Michigan to simulate the coupled 



Sun-Earth system ( Toth et al. 2004, 2005). One of the main modules within the SWMF 
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is the BATS-R-US MHD code. The SWMF is a structured collection of software building 
blocks to develop Sun-Earth system modeling components, in order to couple and assemble 
them into applications. The framework was designed to have "plug-and-play" capabilities, 
and presently it links together nine inter-operating models of physics domains, ranging 
from the surface of the Sun to the upper atmosphere of the Earth. Tying these models 
together gives a self-consistent whole in which each region is described by a world-class 
model, and those models communicate data with each other. The SWMF has achieved 
faster than real-time performance on massively parallel computers, such as the NASA's 
Columbia supercluster. 

4. Simulation Results for CR2107 and Comparison with Observations 

At the rising phase of solar cycle 24, the Sun is becoming more and more active. In 
this study, we focus on the Carrington Rotation 2107 (From 2011 February 16 to 2011 
March 16). During this carrington rotation, a series of interesting events happened on 
March 7 in NOAA 11164. First, an M3.7 flare occurred around 20:00 UT, followed by 
a very fast CME with speed ~2200 km s~^. A gradual SEP event was observed at 1 
AU, which suggests that the CME-driven shock may play an important role to accelerate 
particles. Also, gamma-ray emission above 100 MeV was detected by Fermi/LAT and 
lasted for ~8 hr. This is very unusual since the hours-long gamma ray emission from the 
Sun was observed only three times in the past. Therefore, the flare-related acceleration is 
extremely strong in this event. By simulating this event and validating the results with 
the observations, we can achieve a better understanding of the physics behind the varies 
phenomenon from the Sun to 1 AU. 
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4.1. Steady State Solar Wind and Validation 

To simulate the CR2107, a data-driven boundary condition is used for the inner 
boundary magnetic field. We use the SDO/HMl synoptic magnetogram in this study 
(see Fig. 2). Since the uncertainty of the polar region field, a correction is made to the 



polar field {Sun et al. 2011). To extrapolate the initial potential field, Finite Difference 



Iterative Potential-field Solver (FDIPS) is used ( Toth et al. 2011). The resolution of the 
magnetogram is 3600x1440. 

A spherical grid is used in the simulation. The simulation domain reaches 20 R©. In 
total, 1.3x10^ block is used with 1.2x10^ cells. Adaptive mesh refinement is performed to 
resolve the heliosphere current sheet (HCS). The smallest cell reaches ~10~^ R0, while the 
largest is ~0.8 R©. The radial velocity at X=0 plane is shown in Fig. 3. The fast and slow 
winds are evident in that plot. The fast wind is ~600 km s^^. The slow wind is ~300 km 
s-i. 

In Figs. 4-5, zoom-in velocity and plasma beta with magnetic field lines are shown 
near the Sun. The velocity pattern shows more complex structures comparing with 
the solar minimum condition. The magnetic field show many structures (e.g., streams, 
pseudo-streams). From the plot of the plasma beta, we can see the different regions of the 
corona have very different plasma beta. The polar regions with open field lines have small 
beta ~0.01 near the Sun and increase outward. The high beta regions show the location of 
the current sheets in which the magnetic field is very small. 

In order to evaluate the new steady state solar wind model near the Sun, we compare 
the model output to the electron temperature and density derived from EUV images of 



the Sun by using the DEMT method ( Vdsquez et al. 2010). In general, the DEMT method 



uses a time series of EUV images under the assumption of no time variation and uniform 
solar rotation to derive 3D emissivity distribution in each EUV band. By Local Differential 
Emission Measure (LDEM) analysis, the 3D distribution of the electron density and 
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temperature can be obtained. The DEMT method assumes the plasma is optically thin. In 
this study, we use three bands of EUV observation (171, 193, and 284 A) from SDO/AIA. 

In Figs. 6-7, we show the ratio between the model and DEMT output for the electron 
density and temperature. The sphere is at r=l.O5R0. The ring is between r=l.O35R0 and 
r=1.225R0. For the density, we see a good agreement between the model and the DEMT 
output for most regions. The active regions and polar regions in the model have relatively 
smaller density (Note that the DEMT method has larger uncertainty for the active and 
polar regions). For temperature, we find the model and DEMT have better agreement for 
higher altitude. The model temperature at r=l.O5R0 is lower than the DEMT by a factor 
of ~2. 

In Fig. 8 we provide a direct comparison of the observed EUV images with those 
synthesized from simulations. For better visibility we marked the active regions and 
coronal holes. One can see that our numerical simulations well reproduce the observed 
morphological structures, thus confirming the physical reasonings of the new global coronal 
model. 

5. Conclusion and future work 

At the present stage the global coronal model which is the heart of the solar-heliosphere 
model in the SWMF does not employ any longer the ad hoc heating function, with no 
significant loss in the results quality. For the Carrington rotation CR2107 the simulation 
results are compared with observations. The uncertainty ranges in the model are 
comparatively narrow, except for the reflection coefficient and the uncertainties may be 
further reduced in the course of more thorough validation with observations. The contrast 
ratio in the synthetic images will also be improved. 

Another direction for our further investigation is to improve the solar wind model in 
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the way described in sub-subsection 2.4.2. The reflection coefficient at the time is the most 
uncertain parameter and the only point which is certain about the reahstic reflection is 
that it should account for the wave interaction with reahstic gradients of the magnetic field 
and the coronal density. The solar wind propagation to lAU and the results comparison 
with the in situ observations will be presented in the accompanying publications. 
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Fig. 1. — Typical scales of the transition region: the heat conduction scale (blue), Lt, and 
the gravitational height (red), Lg. 
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Fig. 2. — The synoptic magnetogram of Carrington Rotation 2107 from SDO/HMI. The 
saturation value of the magnetic field is set to be 200 Gauss in order to show the active 
regions more clearly. 



- 34 - 




Y[R] 



Fig. 3. — The radial velocity of the steady state solar wind at X=0 plane. The black boxes 
in the Figure show the blocks and the white boxes show the cells. The thick black line near 
the equator shows the location of the heliospheric current sheet. 
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Fig. 6. — Comparison between the SC and DEMT output near the Sun for electron number 
density. The inner ring shows the ratio between the model and DEMT output from 1.035 
Rq to 1.225 Rq and outside background shows the model output. The iso-surface of the 
Sun shows the same ratio at R = 1.05 -R0. 




Fig. 7. — Comparison between the SC and DEMT output near the Sun for electron tem- 
perature. The inner ring shows the ratio between the model and DEMT output from 1.035 
Rq to 1.225 Rq and outside background shows the model output. The iso-surface of the 
Sun shows the same ratio at R = 1.05 Rq. 
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Fig. 8. — The comparison between observations and synthesized EUV images of the modeL 
Top panels (left to right): Observational images from SDO AIA 131 A, SDO AIA 335 A, and 
STEREO A EUVI 284 A. The observation time is 2011 March 7 20:00 UT. Bottom panels: 
synthesized EUV images of the model. The active regions and coronal hole are marked both 
in the observational and synthesized images, to demonstrate the good reproducibility of the 
observed morphological structures in our simulations. 
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